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by 

T0  To  Williams,  Leo  Robinson  and  7.  L„  Hanson 

n 

t  INTRODUCTION 

The  purpose  of  this  interim  report  is  to  present  prel iminarv  analyses  of 
the  hydrologic  data  col  looted  to  date  on  Lon©  Nan  Coulee.  The  results 
obtained  are  not  considered  to  be  at  all  definitive,  but  onlv  show  trends  which 
can  b©  interpreted  qualitatively.  The  current  status  of  the  investigation 
and  the  course  of  action  by  which  concrete  quantitative  resul ts  may  be 
achieved  are  discussed. 

The  Drainage  Correlation  Research  Project  is  an  investigation  of  peak  discharges 
front’  small  Montana  watersheds*  The  study,,  which  is  being  conducted  by  the, Department 
of  Civil  Engineering  and  Engineering  Mechanics  at  Montana  State  University,  Is  spon¬ 
sored  by  the  Montana  State  Highway  Department  and  the  Bureau  of  Public  Roads*  Four 
small  watersheds  in  Montana  were  selected  In  June  1963  for  intensive  hydrologic  in~ 
stru mentation  and  study*  Precipitation  and  stream flow  data  have  been  coil  acted  from 
the  four  ivatersheds  since  August,  1963*  Additional  instrumentation  Cwind  speed  and 
direction,  soil  and  air  temperature,  and  soil  moisture)  was  added  in  July,  1964* 

Table  I  has  a  fist  of  the  types  and  numbers  of  instrumentation  on  each  of  the  four 
watersheds « 

The  watersheds  are  described  in  detail  in  Interim  Report  0i  9  "Selection  of 
Small  Watersheds/'  May  1965*  Details  on  the  instrumentation  are  given  In  Interim 
Report  &2S  "Precipitation  and  Strearaffow  Instrumentation  for  Smal I  Watersheds,"  and 
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Interim  Report  $3 ,  ^Weather  Stations  for  Small  Watersheds/'  July.  1965. 

Since  the  inception  of*  this  study,  several  runoff  events  of  considerable  signifi¬ 
cance*  including  the  unusual  flood  of  June  1 964-,  have  occurred  at  Lone  Man  Coulee.  A 
number  of  hydrologic  techniques  have  been  applied  to  the  data  collected  thus  far  on 
the  watershed,  and  to  data  collected  by  the  U.S.  Weather  Bureau  and  U.S.  Geological 
Survey .  This  report  describes  tho  current  status  of  the  investigation  with  regard  to 
Lone  Man  Coules. 

II  DESCRIPTION  Of  WATERSHED 

Lon©  Man  Coulee  watershed  is  located  five  miles  south  of  Velier  and  about  two 
miles  south  of  Lake  Francis  in  Pondera  County .  It  is  a  long,  narrow  basin  C6  miles 
long  by  2*5  miles  wide)  covering  14.1  square  miles.  The  basin  is  an  area  of  rolling 
hills  with  rock  outcrops  over  some  of  the  steeper  hills.  Most  of  the  hillsides, 
however,  have  fairly  Sight  slopes. 

There  ar©  two  main  streams  which  drain  the  basin.  Miller  Coulee  drains  about 
two  square  miles,  with  Lone  Man  Coulee  draining  the  remainder  of  the  watershed.  There 
are  two  or  three  small  stock  ponds  on  part  of  the  Miller  Coulee  drainage  area. 

The  length  of  the  main  stream  CLon©  Man  Coulee}  is  7.21  miles.  The  length 
along  the  min  stream  to  the  point  nearest  the  mass  center  of  the  watershed  is 
4.55  miles  CLcaK  The  mean  slope  of  the  main  stream  is  0.6S  percent. 

Watershed  elevations  rang©  from  3775  to  4185  feet,  with  a  mean  elevation  of 
3910  feet. 

Approximately  75  percent  of  the  area  is  in  strip-wheat  farming.  The  steeper 
slopes  are  covered  with  native  grass. 

The  soils  in  the  area  are  predominantly  loams.  A  generalized  brQ&k~$G*m  iss 
50  percent  Morton  loam,  25  percent  Fairfield  loam,  and  25  percent  Scobev  loam. 

Ill  HYDROLOGIC  DATA 

Precipitation  data  is  obtained  from  standard  recording  and  non -recording  rain 
gages  v/hich  meet  Weather  Bureau  specifications.  The  recording  gages  are  dual  traverse, 
12-inch  capacity  weighing  type  gages.  Other  meteorological  information  is  recorded 
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on  strip  charts*  as  described  in  Interim  Report  $3. 

A  continuous  record  of  water  stage  has  been  obtained  since  installation  of  a 
water  stage  recorder  at  the  mouth  of  ths  watershed  in  August  8963.  A  stage-discharge 
relationship  has  been  developed*  to  convert  stage  readings  into  discharges*,  The 

stage-d ischar go  curve  is  shown  in  Appendix  A. 

AH  data  Is  reduced  to  digital  form  and  recorded  on  special  sheets*  from  which 
the  values  may  fee  later  transcribed  to  computer  cards .  Reduction  of  data  from  strip 
charts  to  the  sheets  has  been  kept  up  to  date.  Rainfaf 8  histograms*  storm  hydro¬ 
graphs*  and  unit  hydrographs  have  been  developed  from  some  of  the  data. 

Graphs  showing  doily  precipitation  at  each  of  the  four  pro {set  rain- 
gages  are  shown  £n  Appendix  B.  Appendix  B  also  shows  s  tabulation  of  all 
the  peak  flows  which  have  been  recorded  at  Lon©  Man  Coulee. 

IV  FREQUENCY  ANALYSIS 

Ordinarily*  frequency  studies  of  watersheds  are  based  on  ths  analysis  of  rela¬ 
tively  long  streamflow*  or  crest  stage  gage  records.  In  the  present  study*  however* 
records  of  flood  peaks  exist  for  too  few  years  to  make  this  type  of  analysis  feasible. 
Two  alternative  approaches  are  used  below:  Cl)  comparison  with  other  similar  water¬ 
sheds  with  long  records;  and  C25  synthesis  of  frequency  curves  from  the  relation** 
ship  of  rainfall  to  runoff. 

C  8 3  Comparison  w|th  Other  Watersheds 

Comparison  with  other  watersheds  can  be  applied  only  to  a  very  limited  extent 
because  few  watersheds  of  less  than  100  square  mites  area  have  long  enough  stream* 
flow  records  to  be  of  value,  Revert he! ess*  two  near-by  watersheds*  TW©  Medicine 
Creek  and  Badger  Creek*  each  having  30  years  of  record*  have  been  compared  to  Lon© 

Man  Coulee.  The  areas  of  these  watersheds  are  31?  and  133  square  miles  respectively-. 
Figure  I  shows  a  comparison  of  data  front  Badger  Creek*  Two  Medicine  Creek*  and  Lons 
Man  Coulee  on  the  basis  of  discharge  rate  per  unit  of  area.  This  shows  these  water- 
sheds  to  be  quite  comparable  except  for  a  few  extreme  values.  The  U.  5.  Geological 
Survey  has  developed  a  tentative  composite  flood  frequency  curve  for  eastern  Montana 
in  conjunction  with  the  Crest  Stage  Gaging  Program  ([U.S.G.S.  "Interim  Report  on  the 
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FIGURE  1:  Recurrence  Intervals  for  Badger  Creek,  Two  Medicine 
Creek,  and  Lone  Man  Coulee  as  a  function  of  annual 
peak  discharge  per  unit  area. 
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Frequency  and  Magnitude  of  Flloods  in  Eastern  Montana,  "November  19631  which  uses 
ratio  of  the  mean  annual  flood  Ctho  mean  of  the  annual  peak  discharges  of  record^ 
as  th©  ordinate,  and  time  as  the  abscissa.  This  curve  was  not  developed  for 
watersheds  as  fa*"  west  as  Lons  Man  Coulee.  Tho  data  of  Figure  I  have  been  replotted 
in  terms  of  the  ratio  to  the  n® an  annual  flood  and  compared  with  the  U.S.G.S. 
composite  curve  in  Figure  2.  The  frequency  data  for  Badger  Creek  and  Two  Medicine 
Creek  compare  favorably  with  the  U.S.G.S.  curve.  There  are  too  few  data  points 
for  Lone  Man  Coulee  to  determine  whether  any  of  these  curves  are  applicable  to  it. 

( 2)  Synthesis  from  Reinfal I  Frequency  Curves 

Synthesis  of  flood  frequency  curves  from  rainfall  frequency  data  is  another 
possible  approach.  Undoubtedly  the  relationship  between  rainfal I  and  runoff  is 
very  complicated.  However,  for  events  of  long,  hard  rainfall  which  produce  floods 
of  long  return  period,  many  of  the  complicating  factors  such  as  surface  storage, 
non-uniform  rainfall,  infiltration,  etc.  become  approximately  constant  or  negligible. 
For  this  case  a  direct  relationship  should  ©xist  betweeen  the  frequency  of  the 
rain  producing  th©  flood  and  th©  flood  frequency.  A  preliminary  method  to  form  this 
relationship  has  been  developed  by  Mr,  Robinson  and  is  outlined  briefly  below. 

The  widely  accepted  procedure  developed  by  E„  J.  Gumbel  for  applying  extreme 
value  statistics  to  either  precipitation  or  runoff  frequencies  uses  a  record  of 
largest  recorded  value  during  each  year  of  record  of  either 

discharge  rat©  or  precipitation  rate.  In  th©  case  of  precipitation  rate  the  peak 
mean  rate  over  a  short  time  period,  such  as  one  hour,  is  usually  used.  The  mean 
and  standard  deviation  of  the  data  are  formed,  from  which  the  Gumbel  frequency 
line  is  determined.  For  an  outline  of  the  Gumbel  procedure  the  reader  is  referred 
to  a  standard  hydrology  test  (such  as  Linsley,  Kohler  and  Paulhys,  "Applied  Hydro- 
I ogy" ) . 

The  old  rational  runoff  theory  indicates  that  for  long  storms  with  uniform  pre¬ 
cipitation  rate  the  peak  discharge  rate  is  proportional  to  the  precipitation  rate. 

This  suggests  a  proportional  relationship  between  the  Gumbel  fin©  for  rainfall  and 
the  Gumbel  line  for  runoff.  To  estimate  the  constant  of  proportionality,  observe  that 
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RATIO  TO  THE  M£AH  ANf/OAL  FLOOV 


FIGURE  2:  Recurrence  intervals  for  Badger  Creek,  Two  Medicine 
Creek,  and  Lone  Ma,n  Coulee  as  a  function  of  the 
ratio  to  the  mean  annual  flood. 
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from  statistics,  the  sainpie  mean  estimates  the  population  mean  with  easily  deter¬ 
minable  reliability*  The  mean  of  the  annual  peak  discharges  for  a  short  record, 

v' 

perhaps  as  short  as  five  or  ten  years,  will  give  a  good  estimate  of  the  mean  peak 
discharge  for  an  infinitely  long  record.  The  constant  of  proportional ity  can  then 
be  estimated  by  comparing  this  estimated  mean  peak  discharge  with  the  mean  peak 
precipitation  rate  from  the  rainfall  record;  i„©„, 

K  *  YQ  /  YP  __ 

whm§  K  is  the  constant  of  proportional  ity,  Yq  is  the  mean  peak  discharge,  and 
Yq  is  the  mean  peak  precipitation  rate.  The  same  proportional ity  is  assumed  to 
hold  for  the  standard  deviations,  so  the  standard  deviation  of  the  discharge  peaks 
can  be  determined  from 

Sq  «  K  Sp 

where,  Sq  is  th©  standard  deviation  of  the  discharge  peaks,  and  Sp  is  the  standard 
deviation  of  the  precipitation  rate  peaks.  From  the  thus  determined  values  of 

Yq  and  Sq  the  Gumbo I  line  for  th©  discharg©  peaks  is  plotted. 

Some  questions  concerning  the  val idity  and  adequacy  of  the  above  method 
have  yet  to  b©  answered,  e.g.,  the  assumption  that  the  use  of  24-hour  interval  for 
rainfall  data  is  valid.  However,  th©  method  appears  to  adequately  reproduce  flood 
frequency  curves  generated  from  runoff  data  on  watersheds  where  longer  runoff  records 
are  available. 

Throughout  this  type  of  analysis  the  question  arises  as  to  how  reliable  an 
estimate  of  the  return  period  the  Gumbo!  line  actually  gives.  Confidence  curves 
provide  a  useful  technique  for  doing  this.  Since  methods  for  forming  confidence 

* 

curves  are  not  commonly  given  in  hydrology  tests,  on©  method  is  outlined  here.  The 
theoretical  basis  for  this  method  is  given  in  Appendix  C. 

The  equation  of  th©  abscissa  of  th©  confidence  curves  is  given  by 

Ycc  *  Yt  —  Cr.s.e.)  CSy)  /^n$“n) 

where  Yf  is  the  ordinate  of  the  Gumbo t  line  Cthe  value  of  th©  discharge  or 

precipitation  rate),  f.s.e*  is  the  reduced  standard  error,  Sy  is  the  standard 

-  ! 

deviation  of  th©  date,  (Tn  is  the  standard  deviation  of  the  reduced  extremes,  and 
n  is  the  number  of  data  values. 
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Gumbel  tabulated 
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Sy  and  n  are  know  from  the  statistical  analysis  of  the  data, 
values  of  n  end  ^"n;  these  are  shown  in  Table  21 • 


TABLE  H  -  Means  and  Standard  Deviations  of  Reduced  Extremes  £ a tier  Gy mb© I ) 


n 

fh 

n 

<r* 

n 

8 

0.90 

20 

!  .06 

80 

\  .19 

fiO 

0.95 

30 

Ml 

90 

l .  20 

12 

0.98 

40 

1  0-l4 

100 

l  .28 

14 

8.01 

50 

f  .36 

130 

1.23 

86 

1 .03 

60 

1  J? 

200 

1  .24 

88 

1  *05 

70 

1  .89 

oo 

1.28 

The  r.s*©.  depends  only  on  the  return  period  £the  ordinate  of  the  confidence  curves) 
arsd  is  tabulates  in  Table  HI* 

TABLE  221  -  Reduced  Standard  Error  fcr  Various  Return  Periods  C after  Gumbel) 


Return  Period 

ny> 

6  oOl 
MO 
le50 
2.00 
20SO 
3.00 
4  .00 


Reduced  Sfd. 
ErrorCr.s.e.) 

8.85 

135 

I  oOS 
l  .18 
8  o34 

8  .37 
I!  076 


Return  Period 
CTr> 

5.00 

7.00 

80.00 

20.00 

50.00 

100.00 


Reduced  StdL 
Error Cr.s.e.) 

1.98 

2.03 

2.68 

3.38 

4.32 

5.03 


Thus  for  any  values  of  return  period,  the  corresponding  values  of  Ycc  can  b© 
calculated.  The  confidence  curves  are  then  drawn  through  these  points.  CFor 
example,  so©  Figure  3).  These  confidence  curves  show  the  interval  of  time 
within  which  a  given  flood  may  b©  expected  with,  a  probability  of  2/3,  Confidence 
curves  for  other  probabilities  may  be  obtained  by  multiplying  the  r.s«e.  by 
an  appropriate  coefficient. 

Application  of  the  method  of  synthesizing  runoff  frequencies  from  rainfall 
frequencies  to  Lon©  Man  Coulee  is  shown  by  Figures  3,4 „5, 6,  and  7.  Figure  3, 
Figure  4  and  Figure  5  show  the  Gumbel  lines  for  precipitation  at  several  U.S. 
Weather  Bureau  stations  near  the  watershed  end  the  goodness  of  fit  of  the  data 
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FIGURE  3:  Recurrence  of  the  annual  maximum  day  of  rainfall  at 
Valier,  Montana. 
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FIGURE  5:  Recurrence  of  the  annual  maximum  day  of  rainfall  at 
Conrad,  Montana. 
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to  th©  2/3  confidence  curves.  The  Gumbo!  line  for  Conrad,  which  is  15  miles 
from  the  watershed,  appears  to  have  a  better  fit  than  the  one  for  Valier, 
which  is  5  miles  from  th©  watershed;  however,  th©  ©umbel  lines  for  the  two  stations 
are  almost  identical  when  plotted  dimension! ess I y  on  the  ordinate.  Figure  6 
shows  the  Gumbol  lines  for  the  same  stations,  and  the  weighted  mean  of  these 
linos  plotted  with  the  ratio  of  each  peak  precipitation  value  to  the  moan  peak 
value  on  th©  ordinate.  This  makes  the  lines  directly  comparable.  The  weighted 
mean  I  in©  *m$  formed  from  the  individual  station  i ines  by  weighting  them  by  th© 
distance  from  each  station  to  th©  center  of  th©  watershed.  CSe©  Table  SV> 

The  weighted  mean  tins  was  transferred  by  using  th©  equations  given  above  to 
estimate  Sq  (So©  Table  IV)*  The  val&gs  of  Yq  and  Sq  w@p©  used  t©  obtain  th© 
flood  frequency  I  in©  shown  in  Figur©  7.  The  data  points  for  Man  Gout© 

and  the  U.S.6.S.  composite  curv©  ar©  also  plotted  for  corrparison.  This  shows  that 

/ 

the  transformed  rainfall  curve  predicts  a  smaller  discharge  for  Song  return  periods 
than  do  th©  methods  based  on  other  watersheds. 

There  is  not  a  long  enough  record  available  at  Lone  Nan  Coulee  to  determine 
whether  th©  extreme  value  technique  is  more  reliable  than  a  simple  comparison 
with  other  watersheds.  However  both  techniques  point  to  the  fact  that  th©  flood  of 
June  8964  C  th©  extremely  high  data  point  on  Figur©  I  and  Figur©  7)  was  an 
extremely  unusual  ©vent  with  a  very  tong  return  period. 

V.  MULTIPLE  REGRESSION 

Hu I  tip I©  regression  is  a  statistical  process  for  obtaining  the  best  pre¬ 
diction  of  a  dependent  variable  with  several  independent  variables.  The  pro¬ 
cedure  not  only  determines  coefficients  for  th©  independent  variables,  but  also 
gives  sufficient  information  to  tost  th©  relative  importance  of  each  variable. 

The  square  of  th©  multiple  regression  coefficient  CR2)  also  indicates  th© 
goodness  of  fit.  An  R2  of  1,00  indicates  that  all  experimental  values  ar© 
predicted  ©xactly  by  the  computed  equation.  The  greater  the  deviation  between 
actual  and  predicted  values  the  further  R2  deviates  from  1.00.  Another 
indicator  of  the  goodness  of  fit  is  the  standard  error  of  estimate,  which 
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Table  IV;  Computations  for  weighting  the  rainfall  frequency  lines  and  transforming  the  weighted  I 
to  the  flood  frequency  line* 
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FIGURE  6:  Dimensionless  plot  of  the  recurrence  of  the 
annual  maximum  day  of  rainfall. 
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FIGURE  7:  Flood  frequency  curve  transformed  from  rainfall 
frequency  curves  for  Lone  Man  Coulee  compared 
with  the  U.S.G.S.  composite  curve. 


Table  V  „  Bata  for  multiple  regression  of  snowmelt  ©•vents  in  1965 
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indicates  the  deviation  ot  ms  data  about  the  predicted  vat  us  for  any  constant 
set  of  variables. 

As  with  most  statistical  procedures,  the  larger  the  number  of  observations 
the  greater  the  reliability  of  the  prediction  equation.  Since  there  have  been 
only  8  runoff  events  caused  by  rainfall  on  Lons  Han  Coulee  for  which  alt  of  the 
cf  imata logical  data  are  available,  no  attempt  was  made  to  obtain  a  multiple 
regression  equation  from  the  rainfall -runoff  data.  There  were  14  runoff  events 
caused  by  snowmelt.  Although  this  is  still  a  small  number  of  ©vents  for  obtaining 
a  good  prediction  of  peak  flow,  a  multiple  regression  analysis  was  made. 

A  list  of  the  variables  with  definitions  of  each  is  given  in  Tabl©  V.  These 
variables  were  values  which  were  measured  directly  at  the  weather  stations. 

Tlie  measurable  watershed  characteristics  wore  not  included  as  variables  since  these 
would  b©  constants  when  studying  only  one  watershed.  Any  snowmal t  ©vents  with  a  peak 
flow  equal  to  or  less  than  5  cfs  were  left  out.  As  more  data  becomes  available  this 
limit  might  be  increased  to  eliminate  the  small,  frequent  runoffs  and  thus  emphasize 
the  larger  ©vents. 

A  correlation  matrix  is  also  computed  as  part  of  the  multiple  regression 
program.  Thus  if  two  of  the  variables  are  highly  correlated,  one  can  be  left 
out  to  eliminate  the  resulting  complications  in  the  multiple  regression.  A 
correlation  value  of  0.879  was  obtained  between  the  soil  moisture  reading  at 
3  in.  and  the  soil  temperature.  All  the  rest  of  the  correlations  mro  below 
0.800.  Some  of  the  other  high  values  were: 


%  and  M|8 

0.796 

H3  and  Hg 

0,780 

S  and  Hg 

0.763 

AH  of  the  rest  were  below  0.75.  Sine©  non®  were  above  0.90,  the  independent 
variables  wore  not  highly  correlated. 

An  initial  set  of  regressions  was  run  with  the  discharge  as  a  function 
of  each  individual  variable  and  with  various  combinations  of  the  independent 

variables.  The  highest  values  CO. 527  was  obtained  using  all  9  variables; 
however,  the  stendadd  error  of  estimate  was  122  cfs.  By  eliminating  some  of 
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fhs  variables  a  smaller  standard  error  ms  obtained  with  a  s Sight  decrease  in 
the  ccvrelaticn  coefficient*,  The  equation  of  best  fit  is 
Q  *  -163 ,43  *  5„37  h\2  -  3*99  Mg  *  6.08  V 

The  volue  is  0*48  and  the  standard  error  of  estimate  is  90  cfs* 

These  are  stilA  unsatisfactory  values  for  R2  and  the  standard  error;  but  as  more 
data  points  beccjv>  available  a  better  correlation  should  be  obtained *  Other 
variables  may  be  ati<tad«  An  analysis  of  available  snow  surveys  would  give  an 
approximation  of  water  mailable  for  runoff.  A  change  in  temperature  or  a  ratio 
of  A{2  end  A24  might  be  mxv  appropriate  than  a  temperature  average,  or  degree-days 
may  be  the  most  significant  temperature  factor,.  Percentage  of  snow  cover 
might  also  be  a  factor  to  bo  tried. 

Several  multiple  regressions  wera  frieC  using  the  logarithms  of  the  data* 

The  results  war©  not  as  good  es  the  equation  shown  above » 

VI .  RUNOFF  EVENTS 

C 1 I  Distribution  of  Events 

The  runoff  ©vents  have  been  divided  Ini©  rainfall  events  and  snowmelt  events * 
An  event  is  defined  as  the  recording  of  a  crest  stage  peak  or  a  definite  rise 
and  fall  of  the  continuous  record  trace*  Rainfall  ©vents  also  include  snowmelt 
aecastfpanfed  by  rainfall,  since  the  division  was  made  only  as  to  the  type  of 
precipitation  b&sed  on  temperature  and  firm  of  year* 

Hi©  breakdown  of  runoff  events  by  months  is  shown  in  Table  VI*  Of  th©  total 
number  of  events  46$  were  snowmelt  ©vents  end  54$  were  rainfall  events* 

*2)  M.stHbi^tloo,_of  Snow  Cover  and  land.  Use, . 

Aerial  photographs  of  the  watershed  have  been  obtained  for  evaluating  snow 
cover  and  land  use*  An  analysis  of  the  three  sets  of  pho*bs available  have  made 
it  possible  to  determine  the  percent  of  snow  cover  and  the  percent  of  the  area 
in  various  land  uses*  Those  data  are  shown  in  Tables  VII  and  VIII* 
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TABLE  VU  -  l AMD  USE  BY  PERCENT 


Date 

Stubble 

Suramar  Fal low 

Ranqe  Lend 

Other  Cropland 

1-28-64 

42.3 

43.0 

8.6 

6.2 

2-1  2-65 

42.2 

43.0 

8.6 

6.2 

2-23-66 

42.4 

41  .7 

8.7 

1,2 

TABLE  VIII-  PERCENT  OF  LAND  IN  VARIOUS  USES  COVERED  WITH  SHOW 

Summer  Range  Other 

Date  Stubble  Fallow  Land  Crop  Land  Total 


1 -28-64 

97.3 

18.9 

66.6 

84.9 

60.1 

2-12-65 

100.0 

100.0 

100.0 

100.0 

100.0 

2-23-66 

83.2 

65.0 

63.6 

69.1 

72.7 

VII  OMIT  HYDROGRAPH 

On©  of  the  nsost  frequently  used  techniques  in  the  prediction  of  flood 
magnitudes  is  the  unit  hydrograph.  The  principle  of  the  unit  hydrograph, 
which  was  first  ©^pressed  by  L.  K.  Sherman  in  1932,  is  based  on  the  hy¬ 
pothesis  that  identical  stones,  with  the  same  antecedent  conditions,  produce 
identical  hydrographs.  Under  the  unit  graph  theory  all  hydrographs  resulting 
from  rainf&l Is  of  a  given  duration  have  the  seme  tin®  base.  If  the  rain¬ 
fall  distribution  in  the  storms  is  similar  with  respect  to  time  and  area, 
the  ordinates  of  each  hydrograph  will  b©  proportional  to  its  respective  vo'uma 
of  runoff.  The  unit  graph  is  the  hydrograph  resulting  from  a  runoff  of  1.00 
inch  over  the  entire  watershed.  One©  a  unit  graph  has  been  ps&rpared,  a  storm 
of  any  size  can  b©  considered  and  a  hypothetical  hydrograph  obtained,  so  long 
as  th©  storm  is  of  the  same  duration  as  that  for  which  the  unit  hydrograph 
was  prepared. 

Four  of  the  rain  storm-induced  runoff  events  which  hav©  occurred  at  Lone 
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Man  Coule©  a r©  of  sufficient  magnitude  to  warrant  analysis  by  the  unit  hydro¬ 
graph,,  The  events  for  which  unit  graphs  were  drawn  occurred  on  May  3,  1964, 

June  8,  8964,  June  87,  1965,  and  June  25,  8965.  Hydrographs  and  unit  graphs 
for  these  events  are  shown  in  Figures  8,9,80,  and  II,  These  figures  also 
show  rainfall  8  histograms  obtained  from  recording  rain  gages  on  the  watersheds. 

C I )  Event  of  May  3,  8964 

The  May  3,  1964  event  was  the  result  of  a  rain  storm  of  between  44  end 
45  hours  duration,  the  storm  was  apparently  fairly  evenly  distributed  over  the 
watershed,  as  evidenced  by  the  rainfal I  histograms  from  the  Teren  ranch  (near 
the  mouth  of  the  watershed)  and  the  Geiger  ranch  Cat  the  headwaters  of  the 
watershed). 

Total  rainfall  on  the  watershed  was  5.08  inches  (arithmetic  average  from 
the  4  raingages  on  the  watershed)  or  5.82  inches  (using  th©  Th lessen  method  of 
weighting  the  stations).  Since  the  hydrograph  shows  two  distinct  peaks,  th© 
possibility  of  considering  th©  precipitation  as  two  separate  storms ,  of  28 
and  8  6  hours  duration  was  considered.  About  0.40  inch  fell  sporadically  before 
and  after  th©  main  storms .  3.85  inches  fell  during  th©  first  28-hour  burst, 
producing  0,17  inch  of  runoff,  while  0.80  inches  of  rain  fel  I  during  th©  16- 
hour  burst  producing  0.23  inch  of  runoff. 

Since  there  was  m  has©  stream  flow  prior  to  the  storm,  n©  separation  is 
necessary  in  construction  of  th©  unit  hydrogrsph.  As  is  shown  in  Figure  8,  the  unit 
hydrograph  peaks  are  590  and  840  cfs. 

Basin  lag  (time  from  center  of  mass  of  rainfall  to  hydf-ograph  peak)  is 

10.4  hours  for  the  1st  peak  and  8  8.2  hours  for  the  2nd  peak.  Basin  lag  in  each 
case  is  less  than  storm  duration  tims. 

Soil  moisture  measurements  are  not  available  for  the  period  preceding 
the  event.  However,  examination  of  antecedent  precipitation  would  lead  to  the 
assumption  that  the  soil  should  have  been  fairly  wet  at  th©  beginning  of  the 
storm.  Except  for  some  mixed  rain  and  snow  ending  Apr!  8  6,  there  had  bean  no 
precipitation  until  ApriS  22,  On  that  date  the  watershed  received  0,59  inches 
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(averag©  of  Toren  and  Geiger  gages);  April  23,  0.32  inch;  April  24,  0.17 
inch;  April  26,  0.03  inch;  April  29,  0,11  inch;  April  30,  0.47  inch,  for  a 
total  of  1 .69  inches  in  the  9  days  prior  to  f4ay  I .  None  of  this  April  rain 
produced  any  runoff. 

(2)  Event  of  June  8,  1 964 

Rain  started  falling  on  Lon®  Wan  Coulee  watershed  at  2  p.ro.  on  June  7. 

CA  shower  earlier  in  the  day  had  amounted  to  less  that  0.10  inch).  In  the  30- 
hour  storm  that  followed,  precipitation  was  5.30  inches  of  rain  {arithmetic 
average  from  4  rain  gages);  5.20  inches  (using  the  Thiessen  method);  or  5.33 
inches  Cby  constructing  isohyets).  Although  intensity  varied  from  0.01  inch 
per  hour  to  0.61  inch  per  hour,  in  general  it  was  fairly  steady,  averaging  about 
0.20  inch  per  hour*  Mr.  R,  A.  Oightman,  U.5.  Weather  Bureau,  has  compared  rain¬ 
fall  intensities  from  several  recording  stations  (including  the  Geiger  station 
on  Lone  Wan  Coulee)  located  over  100  miles  apart,  all  of  which  received  com¬ 
parable  total  amounts  of  precipitation.  He  found  remarkable  agreement,  with 
the  several  stations  each  recording  nearly  identical  rainfall  patterns,  and 
fairly  uniform  intensities  of  0.20  inch  per  hour. 

The  hydrograph  for  this  event  (which  appears  in  Figure  9),  shows  a  total 
runoff  of  I .67  inches.  Slope  area  measurements  mad©  after  the  flood  showed  peak 
flows  of  1460  cfs  for  Lone  Man  Coulee  and  280  cfs  for  Miller  Coulee,  (both  de¬ 
terminations  made  ebov©  their  confluence  a  few  hundred  feet  upstream  from  the 
water  stage  recorder).  There  is  no  wav  of  knowing  whether  the  two  streams  peaked 
at  the  same  time.  The  water  stage  recorder  showed  2  lesser  peaks,  at  I!  a.m. 
and  I  p.m.;  and  the  maximum  stage  at  5  p.m.  on  June  8.  Mr.  Williams  visited 
the  recorder  site  ©t  3:30  p.m.  on  June  8,  and  found  no  evidence  at  that  time 
that  either  stream  ©as  falling.  The  5  p.m.  peak  registered  by  the  recorder  cor¬ 
responds  to  the  time  of  failure  of  the  72-inch  culvert  immediatel y  upstream  from 
the  recorder.  A  considerable  amount  of  water  had  been  ponded  by  the  culvert, 
and  the  peak  discharge  of  2685  cfs  determined  by  slope  area-measurements  down- 
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stream  from  the  recorder  reflects  the  sudden  release  of  the  ponded  water. 

From  the  evidence  available*  it  seems  most  likely  that  the  two  streams  peaked 
simultaneously  about  5  p.m.,  and  that  the  peak  discharge  was  8460  t  280  or 

1740  cfs  o  ' 

*  ■'< 

The  unit  hydrograph  constructed  for  this  event  has  a  peak  discharge  of 

1040  cfs. 

Basin  lag  is  II  hours. 

Antecedent  precipitation  studies  show  that  there  was  no  rain  between 
May  II  and  May  27*  On  May  28  there  was  0.52  inch  and  on  May  29  there  was  0.15 
inch.  Mo  further  rain  was  reported  until  June  7.  It  might  be  concluded  that 
there  was  some  soil  moisture  deficiency  at  the  start  of  th©  storm  on  June  7. 

'  O)  Evant  g±  iurre  J2»  1965 

A  series  of  rain  storms  on  June  15-17*  1.965  produced  a  peak  of  110  cfs  on 
June  17.  Although  the  intensity  varied  greatly  with  time*  the  rain  was  very 
evenly  distributed  over  the  watershed.  Th©  first  7 -hour  burst  left  0.88  inch 
at  the  Geiger  station*  and  0.89  inch  at  the  Toren  station.  Fight  hours  after 
the  first  storm*  a  second  1 1 -hour  storm  began,  this  time  recording  1.34  inches 
at  th©  6sig©r  station  and  1 .38  inches  at  Toren^s.  This  storm  was  fol  lowed  by 
a  2-hour  break,  after  which  a  3rd  storm  left  0.22  inch  at  Geiger's  and  0.40 
inch  at  Toren's. 

Total  rain fal I  for  th©  35  hour  period  was  2.45  inches  (arithmetic  average 
from  4  raingages >0  or  2*42  inches  (Thiessen  method).  Omitting  th©  first  7-hour 
stOrro,  the  precipitation  was  I  =67  inches  (average  of  Toren  and  Geiger  stations) 

Separation  of  hydrographs  for  th©  individual  bursts  has  not  bee?)  attempted 
The  hydrograph  for  these  storms  (Figure  10)  indicates  a  total  runoff  of  0.17 
inch.  The  unit  hydrograph  has  a  peak  of  650  cfs.  Base  flow  of  3  cfs  preceding 
the  storm  was  subtracted  out. 

Basin  lag,  if  tine  first  seven-hour  storm  is  neglected,  is  81.7  hours. 

A  study  of  antecedent  precipitation  showed  that  th©  watershed  received  fro 
quent  showers  during  th©  first  week  in  June!  there  was  no  further  precipitation 
until  June  12  when  0.46  inch  fell  i average  of  Geiger  and  Toren  gages).  There 


OOOT 


1^4,  S' 


was  0.10  inch  Hune  14.  Soil  moisture  readings  are  available  for  the  period 
preceding  this  event,  and  are  discussed  below. 

£4)  Event  of  June  25-26.  1965 

Rain  falling  in  several  short  bursts  produced  two  peak  discharges,  on  June 
25  and  June  26,  1965.  The  first  burst  at  the  Geiger  ranch  produced  0.64  inch 
of  rainfall  in  one  hour.  The  total  rainfall,  spread  over  a  49-to  53«hour  period* 
was  1.94  inches  ^average  of  4  stations),  or  1.82  inches  CThiessen  method) . 
Considering  the  1st  storm  as  starting  at  noon  on  the  24th  and  ending  at  2  a.m., 
on  the  25th,  0.69  inch  fell  C average  of  Toren  and  Geioer  stations).  From  2  a.m. 

on  the  55th,  throuoh  6  a.m.  on  the  26th,  0,78  inch  fell , 

The  hydrograph  tor  this  event  CFigure  II)  shows  a  total  runoff  of  0.31 • 
inches,  with  0.16  inch  in  the  first  rise,  and  0.15  inch  is  the  second.  Peak 
discharges  of  200  cfs  and  65  cfs  were  recorded.  These  correspond  to  unit  hydro¬ 
graph  peaks  of  1240  and  425  cfs. 

In  determining  the  basin  lag,  precipitation  falling  between  noon  on  June 
24  and  2  a.m.  on  June  25  Ctime  of  the  first  peak)  was  considered  as  the  first 
storm.  Basin  lag  was  9.4  hours.  Precipitation  falling  between  2  a.m.  on 

June  25  and  6  a.m,  on  June  26  Ctimo  of  the  second  peek)  was  considered  as  the  first 

storm.  Basin  lag  was  9.4  hours.  Precipitation  falling  between  2  a.m.  on 

June  25  and  6  a.m.  on  June  26  Ctime  of  th©  second  peak)  was  considered  the 

second  storm.  Basin  lag  was  12.8  hours. 

The  ©vent  of  June  17,  1965  C which  was  discussed  above)  preceded  the  June 

25-26  event  by  only  8  days.  Light  showers  were  reported  on  June  19,  followed 

by  no  precipiation  until  June  23.  The  Geiger  ranch  received  0.22  inch  during 
the  night  of  June  23-24.  The  Field  ranch  reported  0.32  inch  in  the  24  hour 

period  ending  at  noon  on  June  24.  Soil  moisture  conditions  are  discussed  below, 

(5)  sen,  texture  MaaaEasanla 

Soil  moisture  measurements  have  been  made  at  the  Toren  and  Geiger  weather 
stations  since  July  1964,  The  measured  soil  moistures  together  with  th©  re¬ 
corded  rainfall  at  these  two  stations  is  shown  in  Figures  12  end  13  for  the 
runoff  events  of  June  17,  §965  and  June  25^26,  1965. 
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The  soil  moisture  is  measured  with  gvpsura  blocks,  which  require  cal  - 

ibration  to  determine  moisture  percentage.  Th©  cel ibration  curves  for  the 
gypsum  blocks  used  are  not  exact;  hence,  the  values  shown  in  Figures  12  and  13 
may  not  bo  very  accurate,  although  they  ere  comparable  relative  to  on®  another. 
It  is  relatively  safe  to  assume;  however,  that  th©  high  soil  moisture  values 
C20  -  26$  >  fissured  after  Jim©  15  are  near  the  water  holding  capacity  of  the 
soil  . 

Th©  soil  moisture  values  correspond  well  to  the  precipitation,  with  low 
moisture  values  being  recorded  during  th©  first  part  of  June  when  there  was 
little  precipitation.  Th©  soil  during  this  tin®  was  apparently  quit©  dry. 

On  June  16  there  occurred  e  significant  rise  in  th©  moisture  content.  This 
indicates  that  significant  amounts  of  th©  precipitation  were  infiltrated  into 
the  soil  and  did  not  contribute  to  runoff.  Tha  soil  moisture  content  remained 
high  until  the  event  of  June  25-26.  For  this  event  there  was  no  evidence  of 
further  increase  in  soil  rroistur©,  indicating  ©  small  rat©  of  infiltration. 

The  difference  in  the  antecedent  soil  moisture  between  these  two  ©vents  is  re- 

t 

fleeted  in  tha  ratio  of  total  runoff  to  total  precipitation  in  Table  IX,  since 
less  runoff  per  unit  precipitation  occurred  for  the  June  16  ©vent  in  which  more 
infiltration  occurred.  This  soil  moisture  difference  is  also  reflected  in  the 

i 

difference  between  th©  peak  discharge  of  th©  two  events. 

(6)Piscussion  of  tha  Unit  Hvdroqrephs 

As  indicated  ©artier,  unit  graphs  resulting  from  rainfalls  of  th©  same 
duration  are  expected  to  have  the  s m®  time  base  and,  providing  the  rainfall 
distributions  are  similar  with  respect  to  time  and  ere©,  th©  ordinates  of  each 
hydrogreph  should  be  proportional  to  its  respective  vol urn©  of  runoff.  Pertinent 
information  from  th©  ©vents  described  above  is  summarized  in  Table  IX0 
Because  of  the  few  ©vents  available  for  study,  it  is  not  possible  to 
generalize  from  th©  unit  graphs  that  have  been  produced.  Th©  hydrographs  from 
th©  28-hour  storms  hev©  widely  varying  time  bases,  and  the  peak  discharges 
from  the  unit  graphs  sr©  vastly  different.  Less  deviation  is  noted  in  the  two 
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TOLE  IX 


Date 

Duration 
of  storm 

hrs. 

Tim©  Base 
for  hydro¬ 
graph 
hrs. 

Basin  Lag 

hrs. 

Total 

Precip- 

tation 

in. 

Vol  urn© 
Runoff 

in. 

Ratio 
Total 
Runoff 
to  pre- 
cip. 

Peak  discharge,  cfs 
Hydro-  Unit 
graph  graph 

6-25-65 

14 

28 

9.4 

0.69 

0.16 

0.232 

200 

1240 

6-1  7-65 

16 

37 

1  1  .7 

1  .67 

0.17 

0.102 

650 

5-3-64 

16 

33 

1  1  .2 

0.80 

0.23 

0.288 

840 

6-26-65 

28 

53 

12.8 

0.78 

0.15 

0.193 

65 

425 

5-2-64 

28 

29 

10.4 

3.85 

0.17 

0.044 

590 

6-8-64 

28 

36 

1  1  .0 

5.30 

1  .67 

0.315 

1740 

1040 

16  hour  storms.  Only  two  of  the  events  could  be  treated  as  having  resulted  from 

single  storms. 

To  predict  flow  rates  from  storms  of  given  magnitude  using  the  unit 
hydrograph  the  portion  of  th©  precipitation  which  contributes  to  runoff  must 
be  known.  The  part  ©f  the  precipitation  which  does  not  contribute  to  runoff 
includes  interception  and  depression  storage,  evaporation,  and  infiltration 
into  the  soil .  If  these  factors  can  b©  predicted,  the  remaining  portion 
must  b©  runoff. 

With  antecedent  soil  moisture  measurements  mad©  for  only  two  ©vents,  pre¬ 
diction  of  infiltration  is  not  possible  at  this  time.  Interception  and  de¬ 
pression  storage  have  not  yet  been  evaluated,  but  are  expected  to  be  quite  smal I 
in  comparison  with  infil tration.  Evaporation  is  generally  a  small  value  also. 
CONCLUSION 

7fiis  report  has  considered  the  analysis  of  the  data  collected  to  date  on 
the  Drainage  Correlation  Research  Project  at  Lone  S*?an  Coulee.  The  frequency 
analysis  of  the  available  data  shows  pnomis©,  although  only  very  rough  estimates 
of  the  return  period  of  floods  of  given  magnitued  can  be  made.  Th©  multiple 
regression  analysis  of  th©  data  for  snowmalt  ©vents,  showed  only  unsstis factor' 


-  17  ^ 


correlation  between  peak  discharge  rat©  and  the  correlated  variables.  There 
were  not  enough  data  available  to  perform  a  multiple  regression  analysis  for 
the  rainfall  events.  The  unit  hydrograph  analysis  shows  wide  variation  among 
the  tini©  bases  and  peak  discharges  of  unit  hydrographs  derived  from  various 
storms.  The  fact  that  the  basin  lag  remains  relatively  constant  among 
these  storms  reflects  the  basic  soundness  of  the  unit  hydrograph  technique. 

The  main  problem  is  most  likely  the  lack  of  runoff  events  which  closely  fit 
the  assumptions  of  the  unit  hydrograph. 

Each  technique  of  analysis  applied  to  th©  Lone  Han  Coulee  data  shows 
promise  in  so  far  as  it  can  be  applied.  The  only  apparent  problem  is  that  the 
records  obtained  to  date  are  not  long  enough  to  give  results  of  adequate 
rel iabli ity. 
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APPENDIX  B 


RUNOFF  EVENTS  RECORDED  AT  IX>NE  MAN  COULEE 


The  U.  S.  Geological  Survey  installed  a  crest- stage  gage  at  Lone  Man  Coulee 
on  May  27,  1959*  Peak  flows  reported  by  USGS  prior  to  installation,  in  August 
1963,  of  a  continuously  recording  water  stage  recorder  as  part  of  this  investi¬ 
gation,  are  as  follows: 

Discharge  Discharge  per  square  mile 


June  16-18,  1948 

1820  efs 

160  1 

(Slope- area  measurement  prior  to  crest  stage 
gage  installation;  location  of  measurement  not 
clear  **  watershed  indicated  as  11.4  sq.  mi) 

June  26,  1959 

20.2  cfs 

1.43 

June  27,  1959 

18.5 

1.31 

Nov.  21,  1959 

17 

1.21 

Mar.  17,  I960 

23.4 

1.66 

Mar.  1$,  I960 

23.7 

1.68 

Mar.  19,  I960 

23.7 

1.68 

Mar.  20,  i960 

28 

1*98 

Mar.  26,  i960 

10.5 

0.74 

May  1,  I960 

31.5 

2.24 

May  17,  1961 

21.5 

1.52 

Jan.  29,  1962 

15 

1.06 

Jan.  30,  1962 

10.3 

0.73 

Feb.  3,  1963 

20 

1*42 

Mar.  27,  1963 

24 

1.70 

Aug.  19,  1963 

27.2 

1.93 
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JANUARY 
10  15  20  25 


FEBRUARY 
5  lO  15  20  25 


MARCH 
IQ  15  20  25 

r  5 1 J 


APRIL 

lO  15  20  25 


MAY 

lO  15  20  25 


JUNE 

lO  15  20  25 

rt44+r+(-L-!4+-U+L 


JULY 

10  15  20  25 


AUGUST 

5  lO  15  20  25 

r~rTTTrn  rrrrrr 

--L  14 1-144  4-t.  -LZl.i 


SEPTEMBER 
5  lO  15  20  25 
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lO  15  20  25 


NOVEMBER 
5  1015  20  25 
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DECEMBER 
5  lO  15  20  25 
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lO  15  20  25 
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0.65 


lo  15  20  25 
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AUGUST 


5  lO  15  20  25 

SEPTEMBER 
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lO  15  20  25 
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lO  15  20  25 
OCTOBER 
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MAY 
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The  following  tabulation  shows  peal:  flows  recorded  at  Lone  Man  Coulee 
since  installation  of  the  water  stage  recorder: 


March  30,  1964 

45  cfs 

3.19  cfs/sq  in 

April  12,  1964 

18 

1.28 

May  3,  1964 

195 

13*82 

June  8,  1964 

1740 

123.4 

June  17,  1964 

26 

1.84 

Feb,  17,  1965 

390 

27.65 

Feb.  27,  1965 

29 

2.06 

Mar.  31,  1965 

34 

2.41 

Apr.  I,  1965 

68 

4.82 

Apr.  5,  1965 

63 

4.46 

Apr.  6,  1965 

240  ‘ 

17.0 

Apr.  7,  1965 

60 

4.25 

Apr.  8,  1965 

17 

1.20 

Apr.  9,  1965 

ll 

0.78 

Apr.  ll,  1965 

18 

1.28 

Apr.  12,  1965 

18 

1.28 

Apr.  17,  1965 

32 

2.27 

June  17,  1965 

110 

7.80 

June  25-26,  1965 

200 

14.2 

July  12,  1965 

44 

3*12 
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APPENDIX  C 


THEORETICAL  BASIS  FOR  COltFIDMCE  CURVES 


The  reduced  standard  error  (r.s.e.)  for  a  population  of  extreme  values 
(for  large  values  of  sample  size,  n)  is  given  by 


/<x2  N  tr  2  ^  ~  v/ tt 2 / 6  -/-*  /.  13S6  (x  -  y)  it/ Jo  *  U  (k  -  T^2 

where  c<  Is  the  slope  of  the  theoretical  distribution  with  respect  to  the  reduced 
variate  on  its  probability  paper.  V  is  the  population  mean,  o*(x)  is  the 
standard  deviation  of  the  population,  and  x,  is  the  value  of  the  reduced 
variate.  For  the  order  statistic  of  minimum  variance  this  reduces  to  a  function 
of  only  the  return  period  (Gumbel,  "Statistics  of  Extremes,"  pp.  212-216,  227-229). 
These  results  are  given  in  Table  JR . 

The  reduced  standard  error  is  related  to  the  standard  error  of  the  mth 


value  of  y  by 


<r(/j 


pi/ 


(n$.e.)/j7T 


(Gumbel  p.  53).  The  value  of  o<  ,  according  to  the  Gush  el  plotting  procedure  is 


estimated  by 


a  =  <r„  /s„ 


where  sA  is  the  standard  deviation  of  the  reduced  extremes-  and  .G'n  Is  a-  function 
only  of  the  sample  size.  The  equation  of  the  confidence  curves  is  -given  by 

=-•  y-i  -t  <r  ( y„, ) 

where  is  the  value  of  the  variate  on  the  confidence  curve  and  y.  is  the 
cc  u 

value  of  the  variate  on  the  Gusfeel  line. 

Combining  (l),  (2)  and  (3)  gives 

y«  =  )4  *  (’~.s.e.)(sx)/(<r»/w)  7 

which  is  the  equation  given  in  the  text. 

The  above  method  can  be  applied  as  sn  approximation  for  other  plotting 
techniques  besides  the  Gumbel  Procedure,  in  which  case  y,  becomes  the  ordinate 

V 

of  the  plotted  frequency  line. 
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